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Geminga was first detected as a gamma-ray point source by the SAS-2 gamma-ray satellite ob¬ 
servatory and the COS-B X-ray satellite observatory. Subsequent observations have identified 
Geminga as a heavily obscured radio-quiet pulsar associated with a nearby (250 pc) late Sedov 
phase (300,000 year) supernova remnant. The Geminga pulsar is the second brightest source de¬ 
tected by the Large Area Telescope aboard the Fermi gamma-ray satellite (Fermi-LAT) and has 
been frequently advanced as a source of the anomalous excess of cosmic ray positrons reported 
by PAMELA, Fermi-LAT, and AMS-2. It is surrounded by a compact X-ray pulsar wind nebula. 
Observations above 10 TeV by the water Cherenkov observatory Milagro have also revealed a 
diffuse gamma-ray halo around Geminga extending over several square degrees. Since 2007 the 
VERITAS lACT observatory has performed observations of Geminga and the surrounding halo 
region. However, the standard methods of source detection in VERITAS data have insufficient 
sensitivity to angularly extended sources (>0.5°) to reveal a source on the scale of the Milagro de¬ 
tection. In this talk, we describe two approaches being developed to search for angularly extended 
very high energy gamma-ray emission surrounding the Geminga pulsar. 
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1. Introduction 

Geminga is a nearby (250pc) [1] supernova remnant. Due to its proximity to Earth, it may 
play a significant role in the nature of our local interstellar environment. The Geminga pulsar 
has been put forth as the explanation for several astrophysical phenomena including the positron 
excess measured by PAMELA, Eermi-LAT, and AMS-2 [2]. The Geminga supernova itself has 
also been considered as the source of the low density region of the interstellar medium (known 
as the local bubble [3]) surrounding our Sun. Supernova remnants contribute energy to the inter¬ 
stellar medium primarily through the shock front of expanding ejecta and the pulsar wind driven 
by the spin-down luminosity of the neutron star. Pulsar driven winds commonly exhibit VHE 
emission e.g. the Crab pulsar wind nebula (PWN) and the Vela PWN. The VERITAS imaging air 
Cherenkov telescope array has observed the Geminga pulsar location since 2007. Analysis of the 
VERITAS data using standard point source techniques has yielded no significant detection, and 
produced a 99% confidence level limif for a sfeady/poinf source af energies above 300 GeV of 
< 2 * 10^^^photons [4]. A later search for pulsed VHE emission in fhe VERITAS dafa 

also yielded no significanl detection [5]. The Milagro wafer Cherenkov, observafory which has sen- 
sifivify in fhe range 1 fo 100 TeV, has defecfed very high energy emission from a highly exfended 
(abouf 4 °) PWN surrounding fhe pulsar location [6]. The use of esfablished poinf source analysis 
techniques is nof suifable for analyzing data from potentially highly extended sources. With the 
development of techniques capable of analyzing highly extended source data it is possible that we 
may detect extended emission from the Geminga PWN as seen by Milagro. Such techniques will 
also be useful in analyzing other potentially extended sources. 

2. VERITAS 

The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is located at the 
Ered Lawrence Whipple Observatory (ELWO) in southern Arizona (31° 40’N, 110° 57’W, 1.3km 
a.s.L). The array consists of four 12-meter Imaging Atmospheric Cherenkov Telescopes (lACTs) 
and is sensitive to y-rays from 85 GeV to > 35 TeV with an energy resolution of 15-25%. VERITAS 
is able to detect a point source with 1% the flux of fhe Crab Nebula wifhin abouf 25 hours of 
observation and has an angular resolution (68% confainmenf radius) of better fhan 0.1° af 1 TeV. 

When incidenf on fhe Eaith’s afmosphere, y-rays (as well as high energy cosmic rays) induce 
particle cascades, or showers. The charged componenf of fhe shower produces Cherenkov radiafion 
which is observed by fhe individual VERITAS felescope cameras. The resulfing camera images are 
characferized using Hillas paramefers as described in [7]. These images are fhen used fo determine 
sfereoscopic paramefers of individual showers, such as fhe direcfion on fhe sky and core impacf 
position on fhe ground. Cosmic ray showers serve as fhe primary source of background in VER¬ 
ITAS dafa. Selecfion cufs are applied fo fhe combined image paramefers fo drasfically reduce fhe 
number of background cosmic ray evenfs while optimizing sensifivify fo y-ray sources. Addifional 
information regarding fhe VERITAS dafa analysis can be found in [8, 9]. 

3. Analysis strategy 

The ring background model (RBM) and reflecled-region model are widely used fo estimate 
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the level of background emission in lACT data [10]. This background emission typically consists 
of cosmic ray hadrons and electrons which have passed any selection cuts applied to the data and 
are distributed throughout the entire field of view. The RBM method uses a ring around the source 
region to estimate the level of contamination in the region of interest. The reflected-region requires 
observations to be taken offset from the source, then assigns background regions equally offset 
from the camera center as the source region. Figure 1 gives a visual representation of these tech¬ 
niques. The ratio of the area on the sky designated as the source region to the area designated as 
the background region is known as a. A high value of a means that the contribution of the back¬ 
ground estimate’s uncertainty to the final uncertainty in the significance of an observation is larger. 
Because the RBM and RRM require a dedicated background region within the field of view it is not 
possible to obtain suitable values of alpha for significantly extended sources using these methods. 
Furthermore, for significantly extended sources the true morphology is rarely known, leading to 
source photons leaking into the background regions and the obvious drawback of self subtracting 
the source under investigation. VERITAS is currently implementing two analysis techniques which 
have been successfully employed in other VHE experiments to answer the observational challenges 
presented by extended sources. The first of these techniques is the Matched Run Method (also 
known as the ON! OFF method) which selects a second field of view exfernal fo the source and 
uses it for background estimation [11]. The second technique is a 3D maximum likelihood method 
(MLM) similar to the 2D likelihood analysis used in studies of Fennf-LAT data [12]. The MLM 
presented below includes an added dimension based on a gamma/hadron discriminating parameter 
to increase sensitivity to extended sources. 

4. Matched Run Method 

Our first approach to extended source analysis is the matched run method. This method at¬ 
tempts to match a given data observation with a second field of view to be used for background 
estimation. The second field of view needs to have a similar background rate to the ON or data 
observation and is often called the OFF observation. This method has been successfully used by 
previous lACTs such as the Whipple 10m [11]. This approach allows a source filling the entire field 
of view of the instrument to be analyzed. To obtain a well matched field of view, ideally an OFF 
observafion would be taken immediately before or after the data observation is taken, following the 
same path in elevation and azimuth. This is undesirable as it requires twice as much observation 
time compared to the RBM or the RRM. Another way to obtain a good background match is to find 
a different field of view fhat has no evidence of y-ray emission and has a similar background rate. 

A search algorithm has been developed to help identify which fields of view are suifable back¬ 
ground estimation candidates. There are many factors to consider when identifying a possible 
match for a given ON observation, including its elevation angle, azimuthal angle, time of year, 
hardware configuration, night sky background level, weather, etc. Furthermore, the right ascen¬ 
sion and declination of the source can also have a significant effect on the background rate for a 
source, the primary difference being whether it is on or off of the galactic plane. While all of these 
factors effect the background rate of lACTs to some degree, several of them are very significant. 
The elevation angle at which an observation is made has a drastic effect on the number of y-rays 
and cosmic-rays detected, as well as their energy distributions. When searching for a background 
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(b) Reflected Region Method 


Figure 1: The standard background estimation techniques, the Ring Background Method (left) and 
the Reflected Region Method (right). 


match for a given data run, elevation angle is the first parameter considered. The second parameter 
matched in the current algorithm is night sky background. The night sky background level is af¬ 
fected by the right ascension and declination of the observation as well as weather conditions, time 
of year, and elevation angle and azimuth angle of the observation. Closely matching both elevation 
angle and night sky background levels between the source and background observations has been 
effective in finding suitable background regions for several tests thus far. It is important to note 
that we also require the matched run to have the same hardware configuration of the telescopes, 
and prefer to use observations recorded as closely together in time as possible. As tests on other 
sources continue, other restrictions may be required and implemented into the search algorithm. 
At this time no adjustments are being made to the noise levels of the pixels in order match the 
observations more closely, as careful matching of the runs has been sufficient. 

Once a suitable background observation has been found, there are several ways to analyze 
the matched pair. One method is to use one data run and one background run directly. The main 
disadvantage is that only using one background run per ON run results in an a of 1, meaning an 
equal statistical error in the background and signal uncertainty. Using multiple OFF observations 
could decrease the value of a and reduce the statistical uncertainty of the background but comes 
at the cost of either taking additional OFF observations or finding more suitable matched runs. 
Another way to analyze the pair is to use the OFF field of view for a RBM or reflected region 
analysis. Using a ring (as in the RBM) in the OFF field of view can also decrease a and reduce the 
statistical uncertainty of the background estimation. Successful development of this technique will 
be advantageous for several reasons including the ability to analyze previously recorded ON data 
for extended sources, as well as not requiring further OFF observations to be taken in the future. 

The matched run method is under development and has been used in several test cases. Testing 
of the matched run method on fields of view lacking evidence of y-ray emission from standard 
analysis techniques, as well as point source testing, has gone well. Testing on moderately extended 
sources is underway. 
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5. Maximum Likelihood Method 


The maximum likelihood method (MLM) works by modeling the expected distribution of the 
data across a set of observables. This includes modeling both the source and background compo¬ 
nent of the data. Likelihood analyses used in other y-ray detectors, such as the Fenni-GPVT [12] and 
ctools for CTA [13], rely on spatial information alone. For the Fenni-LKV analysis this is sufficient 
as the instrument also includes an anticoincidence detector for reducing background cosmic ray 
events by over 99.97% leading to the dominant triggers being the result of y-ray photons. The re¬ 
sulting contribution to the Fhrm/-LAT background from cosmic ray events is small compared to the 
contribution from background galactic diffuse photons. For lACT detectors, the localization of the 
photons from point sources is sufficient to implement a purely spatial MLM. However, for highly 
extended sources the distribution of y-ray and y-ray-like' events can be quite similar leading to a 
reduced sensitivity of the analysis. To overcome this, the MLM in development for VERITAS data 
analysis includes a third dimension based on a y/hadron discriminating parameter known as mean 
scaled width (MSW). This third dimension is incorporated into the likelihood by including models 
of how both the y-ray and background cosmic ray components of our data are distributed in MSW. 
Details regarding the VERITAS MEM, including information on the likelihood construction, back¬ 
ground spatial model determination, and characterization of the various IREs can be found in [14]. 

MSW is determined as an average of the Hillas width parameters from each telescope image 
normalized by the expected width derived from simulated air-showers with the same integrated 
image intensity and shower core impact distance from the telescope. 


MSW = 



n 


I 


Wj 

< Wi > 


(5.1) 


It can be seen from Eq. 5.1 that the MSW distribution for y-ray events is expected to peak at 
~ 1. Cosmic ray events typically produce broader showers^ resulting in images with much larger 
calculated widths. This leads to a peak at higher values of MSW [7] as compared to y-ray showers 
(see Eigure 2). Of the standard set of selection cuts used to reduce the level of background in 
VERITAS data, a cut on MSW is among the most efficient. Eor a typical VERITAS source analysis, 
events with MSW>L1 are cut. However, in the MEM analysis the MSW range is extended to 1.3 
to provide ample separation between the y-ray and cosmic ray MSW distributions. The distinct 
difference in the shape of the MSW distribution for these two components improves the MEM’s 
ability to distinguish source emission from background as opposed to using a spatial only model. 

The source models are based on the instrument response functions (IREs) relevant to the VER¬ 
ITAS detector. The IREs are combined based on the formulation of Mattox (1996) Eq. 2 [2]. The 
actual construction of the source model is made more tractable by computing the expected spatial 
models in 0.025° x 0.025° bins. Each IRE component is also computed in narrow logarithmic bins 
of true energy^ in steps of 0.05. The final source model for a given spatial bin is fhen compufed 

* y-ray-like events typically consist of cosmic ray electrons and/or hadrons which survive the standard analysis cuts 
chosen to optimize sensitivity to y-ray sources. These events are the main source of background in VERITAS data. 

^This is due to the larger transverse momentum of the secondary pions produced in hadronic cosmic ray showers. 

^true energy is the true energy of an observed y-ray. Observed (or reconstructed) energy is the energy of the photon 
as calculated by the VERITAS analysis. 
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(a) Extended source spatial model 


(b) y-ray MSW distribution 



Mean Scaled Width 


(c) Background spatial model (d) Background MSW distribution 

Figure 2: (a) Extended source spatial model as derived from the Milagro C3 source parameters de¬ 
scribed in [6]. This model is asymmetric due to the observation offset of 1° and the changing sensi¬ 
tivity of the instrument over the field of view, (b) MSW distribution derived from y-ray air-shower 
simulations.(c) Background spatial model and (d) background MSW model derived from VERI¬ 
TAS data. Both background distributions are derived from events greater than 0.4° from known and 
potential y-ray sources. Additionally, the spatial background distribution uses events with MSW < 
1.3 reflecting the range currently used in the MEM analysis. Distributions were derived assuming 
70° elevation, due south, with four telescopes participating in the shower reconstruction. 
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according to 

_ r^tucix itppsr 

S,rciri,j\s) = X Bm,nPmAhMm,n{E'k) / R,nAE,E'k)dE / ^ ’ S{E'\s)dE' (5.2) 

where the components consist of the intrinsic source morphology {B), the instrument point spread 
function (P), effective collection area (A), energy dispersion (P), and source spectrum (S). Spatial 
coordinates are given by f,- E’ represents the true energy of an event, and E the observed energy. 
To account for the smearing of events due to our point spread function, the contriubtion from every 
bin in the map (indexed by m,n) is summed at the position of interest (indexed by ij). Each IRE is 
parameterized for a given set of observing conditions in zenith angle, azimuth, source offset, and 
camera noise level. The IREs are also computed in narrow bins of E’, indexed by k. Integrating 
the computed source model over every bin (iJ) yields an estimate of the total number of observed 
photons in a given dataset. 

To demonstrate how the inclusion of the MSW parameter aids in the detection of extended 
sources, a comparison of the models of an extended source were compared to those of the expected 
background (see Eigure 2). The extended source uses as input the extension and spectrum as 
reported in [6] for the Milagro C3 source. The morphology (B) is assumed to be a Gaussian 
with a=1.3° and the spectrum is modeled as a powerlaw with spectral index of -2.6 and flux of 
37.7x10^^^ (TeV^^cm^^s^^) at 35 TeV. The resulting source model was then derived for the 
VERITAS energy range 0.5 - 1 TeV. When compared to the expected background spatial model, 
the derived extended source model is very similar. However, in the MSW dimension there is a 
significant separation between the two models. 

The presented extrapolation of the Milagro C3 source to lACT energies should only be consid¬ 
ered a test of the MEM. In reality, due to the higher energy threshold and larger uncertainty in the 
event position reconstruction of the Milagro detector, any emission potentially seen at VERITAS 
energies may have a very different morphology and spectrum. 

6. Conclusions 

The Geminga pulsar is an interesting y-ray source with a rich discussion surrounding it in 
connection with our stellar neighborhood. A VERITAS detection of either a point source or an 
extended PWN would contribute significantly to our understanding of the origin and evolution of 
the local interstellar medium. Two analysis techniques are being implemented by VERITAS which 
will aid in the study of spatially extended y-ray emission which may emanate from this region. 
These techniques are both progressing well and should be ready to analyze the VERITAS Geminga 
data in the near future. 
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